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Abstract 

We develop a generalized theory of (meta)equilibrium statistical mechanics in the 
thermodynamic limit valid for both smooth and fractal phase spaces. In the former 
case, our approach leads naturally to Boltzmann-Gibbs standard thermostatistics 
while, in the latter, Tsallis thermostatistics is straightforwardly obtained as the most 
appropriate formalism. We first focus on the microcanonical ensemble stressing the 
importance of the limit t — > oo on the form of the microcanonical measure. Inter- 
estingly, this approach leads to interpret the entropic index q as the box-counting 
dimension of the (microcanonical) phase space when fractality is considered. 
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1 Introduction 



Recent years are witnessing a growing interest on foundational aspects of 
statistical mechanics and thermodynamics in an effort to extend methods 
of standard Boltzmann-Gibbs (BG) thermostatistics well beyond their tra- 
ditional domain of application. Systems exhibiting long-range interactions 
(self-gravitating systems, two dimensional vortices, metallic clusters, etc.), 
fractal memory, spatio-temporal chaos, non exponential relaxation and non 
maxwellian velocity distributions, pose formidable problems to standard BG 
thermostistics or just make it unapplicable. Experimental signatures of these 
systems include, among others, negative specific heats, inequivalence of sta- 
tistical ensembles, phase transitions and self-organization. 

Because of the great success of BG thermostatistics in explaining and pre- 
dicting thermal properties for a wide class of systems, this formalism is a 
cornerstone in which any generalization should be inspired and to which it 
has to be reduced under certain conditions. What are these conditions and 
the most suitable points to introduce generalizations in the standard theory 



have been issues o 
othermo dynamics 



intense debate and controversy. Theories like Hill Nan- 
l[ proceed from standard thermodynamics extending the 
Gibbs-Duhem equation very much in the manner in which Gibbs extended 
the former formalism of thermodynamics giving rise to the extremely impor- 
tant concept of chemical potential. Other generalizations introduce nontrivial 



insights on the standard formalism, as the Hessian of t 



le microcanonical en- 



tropy in studying phase transitions in small systems 21] (rejecting statistical 
ensembles other than the microcanonical out ot the thermodynamic limit). 



Another quite original approach due to Tsallis 



q] gave rise in 1988 to 
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an entirely new field in Statistical Thermodynamics focused on the concept 
of entropy and its appropriate form to deal with complex systems. We have 
shown recently how Tsallis Thermodynamics can be connected to Hill's Nan- 



othermo dynamics 



Tsallis thermostatistics has received wide attention recently, proving quite 
useful in the analysis of nonlinear low dimensional maps and rnultifractals 

besides other many applications including: stellar polytropes 

n 

butions [9], anomalous diffusion 
long-range interactions 
lence 



, ^ ^ „, Levy distri- 
lo| . inelasticity in hadronic reactions jlll ]. 



12| . ionic solutions and fully developed turbu- 
14| . If we consider a continuum spectrum of configurations uj E Q 
{doj = Yli dpidqi) where fl is the phase space volume of the system, Tsallis 
entropy has the following form 

s, = J«f!M!^ (, e ») (1) 

where q is the so-called entropic parameter. It is easy to see that when q tends 
to unity, the well-known Gibbs-Shannon entropy Si = —ki Jj. p{u)ln p{u)du! 
is regained. 

When equiprobability for the F microstates of the phase space is considered, 
Tsallis entropy takes the form 

which is the generalization of the Boltzmann microcanonical entropy (when 
g — > 1, the Boltzmann entropy Si^^^^ = fci InF is regained). 

Maximization of Tsallis entropy, Eq. (1), through appropriate constraints leads 
to q-exponential distributions. These are of great interest to study the power 
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laws which many complex systems exhibit. Tsallis q-exponential distributions 



have been extensively founded through rigorous lines 1^ and have the follow- 
ing form 

cx(l-(l-g)/?E,)T^ (3) 
When q tends to unity the Boltzmann factor is regained. 

The methodology to deal with complex systems is, however, far from being a 
closed debate. Other interesting generalizations have appeared recently based 



in the use of deformed logarithms or kappa distributions 



3- 



m 



17| . Kolmogorov- 



Nagumo averaging [18] and Renyi statistics [19|. A plethora of new entropies 
have also appeared during the last decade, many of them inspired in the Tsallis 
MaxEnt approach 2o| . 

Although Tsallis entropy has many desirable mathematical properties such as 



concavity and Lesche stability |2 



also associated an H-theorem 
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^besides a uniqueness theorem 



22| | having 



23l |. the rationale behind its form and why 
it is so suitable to deal with the thermodynamics of complex systems against 
other entropic forms has not been clarified. Tsallis introduced his entropy 
inspired in the study of multifractals. The entropy appears in a postulate- 
like manner in his first work on the subject 3| and, although its interest 
and beautiful properties (pseudoadditivity, concavity for g > 0, etc.) are soon 
apparent, it would be desirable to have an approach giving some insight in 
the starting point, i.e. the form of the entropy, which is the building block of 
the formalism. 



In this paper such an approach is attempted using concepts in which Tsallis 
himself has been insisting since his above mentioned proposal. The crucial 
issues are: the infinite time limit, nonergodicity and fractality of the phase 
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space. When the latter is thermodynamically important, the phase 



a smooth surface and standard thermostatistics does not apply 



space is not 



2^. Although 



the existence of fractality in the phase space of many-body complex systems 
far from equilibrium is an issue under debate there is plenty of evidence in 

n 

low dimensional systems (see for example ||25i] where the nice example of the 
Galton board is discussed). We show here that fractal analysis of the phase 
space leads quite naturally to Tsallis entropy. The approach allows us to give 



geometry adding insight to previous findings 
to thermal variables and/or their fluctuations. 





m 


I 
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in which q was related 



Although we are concerned also with ensemble theory, our reasoning line here 
is strictly based on the microcanonical ensemble and on how the findings in 
this ensemble translate to the other ones by means of the Legendre transform 
mechanism. In this study Massieu-Planck entropic potentials play an impor- 
tant role in generalizing our results to other ensembles. The interest of these 
quantities at nonequilibrium problems were recognized at the very moment of 
their inception more than one century ago .^J} . These allow for a compact for- 
mulation of generalized statistical mechanics and thermodynamics (see 28| for 
a Jaynes-like formulation of equilibrium standard thermodynamics by using 
these entropic potentials). 



2 Microcanonical measures and ergodicity 

Let us consider a system of particles isolated from its surroundings by rigid, 
impermeable, adiabatic walls. The trajectory of the motion is made in such a 
way that energy is constant at all times (provided that the system is perfectly 
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isolated) . Let us consider a phase space with volume ft containing all possible 
configurations that the system can attain. Each of these configurations is a cell 
of volume ~ h^^ in which the phase space is subdivided. If we want to count 
the total number of cells visited by the actual trajectory of the motion during 
a time T this can be achieved exactly by means of the following mechanical 
quantity 

dt X n (Hiuj) - E) m\m (p. - p^) ^ - (4) 

where Qiit) and Piit) are, respectively, the generalized position and momenta 
coordinates which are the solution of the Q)N Hamilton equations of motion 

r ^ dH 

The phase space visited by the trajectory is thus embedded in the available 
phase space F (the constant energy surface in the volume fl). Note that in 
Eq.(4), the first Dirac delta containing the difference H{p, q)—E is redundant 
if in the actual trajectory specified by (P, Q) energy is conserved. It is written 
here, however, for our convenience, emphasizing that the trajectory lies in 
the constant energy surface. Not all points in the constant energy surface, 
however, are actually points of the trajectory and this is what the other delta 
functions select: only points on the mechanical trajectory do contribute then 
to Eq.(4). 

We now define what we are going to call from here on actual microstates f 
f= lim fdtj Xl'^^8{H{u)-E)\Qm5{p,-P,)5{q,-Q,) (6) 
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We introduce now also the smooth microcanonical measure F or, what we will 
call from here on smooth or Boltzmann microstates 

il[j0'HHiu.)-E) (7) 

We should say that the system is ergodic if F = F. It is clear that the system 
will be ergodic if 

lim r dt^\Qi\\Pi\5{pi-Pi)5{qi-Qi)^l (8) 

T— >oo Jq " 

for each point u — (q, p) e F. This means that each point of the attainable 
phase space (with -^^(q, Tp) — E — constant) becomes, in fact, a locus of the 
trajectory in the long time limit, i.e. the trajectory fills densely the constant 
energy surface. It is important to note that writing here the limit T tending 
to infinity contains a bit of abuse of language since time must be, in any 
case, lower than the Poincare recurrence time to make sure that each phase 
space cell is not visited more than once during the trajectory. The hmit T — > 
oo should be understood as a long time enough to make of a moving point 
in the phase space (a microscopic description of the system) a trajectory 
which is long enough to be considered itself as a macroscopic geometric object 
(if we think in a many particle system at finite temperatures). This is the 
key idea to understand our approach and is applied when the trajectory has 
achieved its maximum spread over the phase space surface. We then do not 
care about of in which phase space point of the trajectory is actually the 
system: the macroscopic trajectory becomes the key geometrical object to 
pass to a statistical description from the microscopic mechanics. We then 
attribute the physical properties of the motion of the system [mechanically, a 
moving point in a high number of dimensions) to an average property of the 
whole trajectory {statistically, a complex geometrical structure embedded in 
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the surface of constant energy). 

Let us see that in the case of an ergodic system the time average (mechan- 
ical average) and the ensemble average (i.e. over the smooth microcanonical 
measure) of a phase function f{P{t), Q(i)) coincide. We define the long time 
average of /, over the trajectory < / > as 

x/(P((), Q(())|4l lAl* (p. - -P.) S (q, - Qi) (9) 

Because of the product of Dirac deltas in the time dependent integrand, we 
have f(P{t), Q{t)) — f{oj), so that / can go out of the time integral. Applying 
Eq. (8) and f = F we obtain 

< / >= f X n m.) -E) = f (10) 

i.e. the time average of the phase function / equals the ensemble average /. 

Ergodicity is the mechanism that allows us to pass from a mechanical descrip- 
tion of the system to a statistical one. We can wonder what should we do if 
the system is not ergodic i.e. if F 7^ F. It is clear, in any case, that the measure 
F is embedded in the smooth measure F composed of all microstates in the 
constant energy surface. If we can still find F as a function of the embedding 
smooth surface F in a form F = A''(F), we can still attempt a relatively simple 
statistical description. In what follows, we will focus in this smooth surface 
considering physical cases in which the function A/'(F) can be constructed. For 
the ergodic case, A/'(F) = F. If the system is not ergodic, we can still have a 
relatively simple statistical description of it if the trajectory is a fractal curve 
embedded in the smooth energy surface of F cells. If the reduced box-counting 
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dimension of such a curve is d {= D/6N so that < c? < 1) and provided that 
;he volume A (oc l/F) of each phase space cell is vanishingly small, we have 



The available phase space is now a fractal object embedded in the Boltzmann 
6A^-dimensional smooth phase space. Clearly, when d = 1 ( 6 N- dimensional 
phase space) we regain ergodicity. 

3 General definition of the microcanonical entropy and particular- 
ization for both smooth and fractal phase spaces 

The first step in passing from mechanics to statistics on the phase space has 
been completed in the previous section. The relationship to thermodynam- 
ics has to emerge now to have thermal variables founded microscopically. In 
doing that we can make the following (necessarily heuristic) reasoning. In a 
nonequilibrium situation, the change in the number of smooth microstates 
(i.e. on the area of the smooth energy surface) can be thought to be directly 
proportional to the number of actual microstates of the trajectory. This can 
be clearly understood if we think on each point of the actual trajectory as a 
source for new smooth microstates. It is in this way that our lack of knowledge 
on the microstate in which the system is can spread, tending to increase. The 
smooth energy surface and the actual trajectory of the system which spreads on 
it are inextricably related to each other. Since we do not know indeed in which 
microstate is the system, we are forced to admit that each microstate behaves 
in nonequilibrium as a source for new ones and these new microstates make 
also the smooth surface to increase. All these considerations can be written 



29 




Af{T) oc r'^ 
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mathematically as follows 



f = .Mr) (12) 



Here a is the rate of smooth microstates production from the actual mechanical 
ones. It is exactly at this point that we introduce thermodynamics by estab- 
lishing the equivalence between this rate and that of entropy production. 

Here S = S/k is the dimensionless entropy. If we replace Eq.(13) in Eq. (12) 
we obtain 

' s-rj^ (14) 



dr AA(r) Ji M{r 

where we have considered in the integration that for Fq = 1 the entropy van- 
ishes. This constitutes a general definition of entropy that, as we are going to 
see, reduces to well known entropic forms after considering models for A/'(r). 
Let us consider first the ergodic case (smooth phase space with A/'(r) = F as 
discussed in the previous section). Prom our definition we obtain the Boltz- 
mann entropy after integration 

5 = lnr (15) 

If, however, a fractal phase space is considered (A/'(r) cx F"^), we have, after in- 
tegration (absorbing the proportionality constant in a generalized Boltzmann 
constant kd so that S — S/ kd) 

Quite remarkably, this form of the microcanonical entropy is similar to that 
of Tsallis, Eq. (2) under the correspondence q = d.li this connection is made, 
our view implies that the entropic parameter q is equal to the reduced box- 
counting dimension d = D/6N of the available fractal phase space. Of course. 
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for Boltzmann systems, D = 6N, d = 1, q = 1 and the Boltzmann entropy is 
regained. These considerations are relevant to many physical systems due to 
the multifractality that thermostatted and low dimensional systems exhibit 
25[ (here we are considering isolated systems in the microcanonical ensemble 
and, therefore, we are focusing in the fractal support of the phase space) . 



This insight allows us to interpret unequivocally the strong coupling regime 



in Ref . 



found in ionic solutions. In [13^ we observed that q tends to vanish 



in this limit. It was previously known that when multivalent ions of the same 



charge are located c 



a Wigner crystal 



ose to a highly charged surface these crystallize forming 



3l|. This means that the available classical phase space col- 



lapses into regions with strikingly lower dimensions. Therefore, D « 6N, 



c? — > and g — > 0, as previously obtained in [13| because of the crystalline 
ordering at the interface. In general, for moderate coupling, there is quasi- 
cristallyne behavior at the interface which is still markedly different to that in 
the bulk (a disordered, uncorrelated phase). The transition between the phase 
space for the bulk and the interface can, hence, be depicted as a complex 
filamentary structure having fractal properties. At low coupling values, the 
liquid in the interface has bulk properties js^, and both, bulk and interfacial 
liquid merge into a single homogeneous and smooth phase. All this variety of 
behavior is explained economically through the index q. The weak coupling 
limit can also be understood analytically following a previous insight in which 
we related q to the interfacial and bulk entropies j6|. 
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4 Generalized canonical ensemble 



The above methodology can be easily extended to the canonical ensemble 
provided that we can characterize the space in which the typical trajectory 
is embedded by a number of configurations J\f{T*) where T* now denotes the 
total number of smooth configurations in the composite phase space surface 
appropiately weighted by the different energies that the different regions of 
the surface have. Following the same principle for the evolution of the number 
of microstates (now weighted) we have 

dr* 



dt 

where 



a*^^{^*) (17) 



.• = t^ (18) 

dt ^ ' 

is the rate of change of the natural entropic potential which in the canonical 
ensemble is minus the Helmholtz free energy F in kT units. This entropic 
potential is related to entropy by means of the Legendre transform mechanism 

F E 

Note that the natural entropic potential correctly accounts for the entropic 
contribution of the phase space cells substracting to them the energy contribu- 
tion which is now a fluctuating variable. If we deflne J- = —F/kT and replace 
Eq.(18) in Eq.(17) we obtain 

dJ^ 1 _ f^} dV 

T 



dv* Af{r*) Ji Af{r*) 

where we have used that for Fq = 1 the entropic potential vanishes. For an 
ergodic system, despite that now regions of the phase space are weighted by 
their differing energy, all them can be attained in the long-time limit so that 
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J\f{r*) = r* and 

:r = inr* (21) 

The total number of "weighted microstates" F* is the usually called partition 
function of the system. For a fractal phase space, although the microstates 
are weighted by energy, the microcanonical restriction coming from a lower 
dimension of the actual phase space still applies and we have after integration 

p*l— d 2 

Our aim now is to calculate the weights in the microstates F* to relate them to 
regions in which the phase space is partitioned. It is important that if subindex 
i denotes regions of constant energy Ei we have 

r* = En (23) 

i 

r = Er^ (24) 

i 

and if for each of these regions Eq.(19) applies 

Ti^Si- PEi (25) 

Now, since the definitions Eqs. (14) and(20) apply also to Si and J^i by re- 
placing these in Eq.(25) we can calculate F* as a funtion of each Fj and the 
energy E^ 

For example, for an ergodic system we have, from the latter equation 

F* = TiC-^^^ (27) 

and, for the nonergodic system with fractal phase space (A''(F*) = F*"*) we 
have 

F* = F, (l - (1 - d)PE,/Tl-'') ^ (28) 
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which coincides with the TsaUis q — exponential form with the only difference 
that /3Ei is here divided by the degeneration of the microstate. Note, however, 
that usually the MaxEnt approach requires the specification of the entropic 
form as well as the appropriate constraints. Here we have not chosen, in prin- 
ciple, any entropic form, nor have privileged any choice for the constraints: we 
have based all our reasoning in the fractal analysis of the phase space and in 
the evolution equation for the number of "smooth" microstates. 

The partition functions for each system can now be calculated from Eq.(23). 
The probability of having the system with energy is given in each case by 
the ratio r*/r*. 

5 Massieu-Planck thermodynamic potentials and generalized ther- 
mostat istics 

In the previous sections we have formulated equilibrium generalized thermo- 
statistics in the canonical and microcanonical ensembles. In general, the mo- 
tion is such that a set of extensive variables Xj are kept constant during the 
trajectory (while fluctuating their conjugate intensive ones yj) and the inten- 
sive variables i/k are also kept constant (and then their extensive conjugate 
ones Xk fluctuate). This always occurs in semi-open systems in which some 
extensive variables fluctuate. We can now define a generalized thermodynamic 
potential X by Legendre transforming the entropy, following a totally analo- 
gous procedure to the sketched above for the canonical ensemble 

X^S-^VkXk (29) 

k 
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This is the thermo dynamical definition of the generahzed potential X. Its 
character of thermodynamic representation is made explicit in its generalized 
differential Gibbs form 



dX^Y.yjdXj-Y.^kdyk (30) 

j k 

from which it is seen that non-environment variables i/j and Xk can be obtained 
from the natural ones by differentiation 

dX dX 

by keeping constant all other variables not involved in the differentiation. In 
Table 1, the sets of variables and the form of X for each statistical ensemble 
are indicated. 

The typical long-time trajectory contains now Af{T) microstates where the 
smooth microstates in the constant X/s surface T obbey the following evolu- 
tion equation when equilibrium is perturbed 

^ = (32) 

which leads to a microscopic definition of X (in a way totally analogous to 
previous sections) as 

'''' 

and we also have from Eq.(29) particularized for a micro- value i of each ex- 
tensive variable Xk (denoted by Xfe,j) 

wk^k jm'^^'^'' ^^^^ 

which allows to know Tj as a function of all the weights and then evaluating the 
generalized partition function from T = Y.i^i (where is to be understood that 
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the number of sums is the same than that of fluctuating extensive variables). 
For an ergodic system, for which A/'(Tj) = Tj, this leads to the generahzed 
Gibbs distribution 

T ^Y.^i^'^"^"^"'' (35) 

i 

and, in the non-ergodic case with fractal phase space (A/'(Tj) = T^) 




6 Conclusions 

The unifying principle of the smooth phase space evolution Eq. (32) allows 
to build a generalized thermostatistic formalism concerning macrovariables 
analogously to that of Hamiltonian mechanics containing microvariables (com- 
parison of Eqs. (5) with (31) makes explicit this analogy: the role played by 
X in macroscopic systems is quite similar to the one played by H in mi- 
croscopic ones). It is to be noted that this principle embodies the maximum 
entropy principle in the microcanonical ensemble as well as the minimum of 
the conventional thermodynamic potentials at equilibrium (a maximum of 
their closely related Massieu-Planck entropic potentials, see Table 1) for each 
ensemble. This is clearly seen in the differential equation which serves as def- 
inition of the entropy Eq. (14), since in the microcanonical ensemble l/A/'(r) 
is directly related to probability. As probability is a positive quantity which 
decreases with increasing the available phase space (of course, this applies far 
from phase separation) entropy is here, therefore, a concave and increasing 
function of its argument. The generalized Massieu-Planck entropic potential 
X describes all the thermodynamics of a given system when known as a func- 
tion of the natural variables and is the key stone, in our view, in connecting 
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the macro to the microworld. 



Our principle for the evolution of the smooth surface of microstates vaguely 
reminds the "differential equations" formulation of q-exponential distributions 



32j |. Here, however, the formulation does not come from mathematical moti- 
vation and/or developing of mathematical identities, but an attempt has been 
made to substantiate it physically. 

We have shown that our picture leads naturally to BG statistical mechanics as 
well as to Tsallis formalism, without need of formulating the latter by intro- 
ducing ad hoc constraints. Tsallis entropy gains advantage in our formulation 
compared to other entropic measures since its physical meaning is here appeal- 
ing and has far reaching consequences beyond any ad hoc formal construction 
with whatever nice given properties. 

We want to acknowledge B. A. Nikolov and J. A. Manzanares for their com- 
ments and critical reading of a previous version of this manuscript. V. G-M. 
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the MCYT (Ministry of Science and Technology of Spain) and FEDER under 
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Table 1 

Environment variables {Xj} (extensive) and {yi} (intensive) and non- environment 
fluctuating ones {X^} (extensive) {yk} (intensive) for each ensemble indicated. {/3 = 
l/kT, TT = p/kT and u = -fi/kT). 
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